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Abstract 

A noncentral x 2 density is log-concave if the degree of freedom is v > 2. We complement 
this known result by showing that, for each < v < 2, there exists X v > such that the % 2 
with v degrees of freedom and noncentrality parameter A has a decreasing density if A < A„ , 
and is bi-modal otherwise. The critical A„ is characterized by an equation involving a ratio 
of modified Bessel functions. When an interior mode exists we derive precise bounds on its 
location. 

Keywords: distribution theory; log-concavity; log-convexity; noncentral distribution; 
special functions; unimodality. 

MSC2010 Classifications: 60E05; 62E15; 33C10. 

1 Introduction 

The central x 2 density with v degrees of freedom is log-concave if v > 2 and log-convex otherwise. 
The noncentral x 2 density with v degrees of freedom and noncentrality parameter A is still log- 
concave if v > 2. It is decreasing if < v < 2 and A is small. When < v < 2 and A is large, 
however, a new shape emerges which some may find surprising: the noncentral x 2 density can 
be bi-modal! 

This note characterizes the noncentral x 2 density as either log-concave, decreasing, or bi- 
modal (the first two categories overlap when v = 2) by delineating the range of parameters for 
each. Our contribution is a mathematical derivation of some facts about the noncentral x 2 that 
are "folklore," i.e., either tacitly assumed or supported by numerical evidence. This happens 
to be an interesting exercise in special functions. Part of the criteria is expressed in terms of a 
modified Bessel function equation, and the method may yield further results (e.g., inequalities) 
concerning modified Bessel functions. 

2 Definitions 

The central x 2 density with v > degrees of freedom is given by 
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The noncentral x 2 density with v degrees of freedom and noncentrality parameter A is defined 
as 

00 e -A/2 ( \\ k 

Pu,x(x) = jp— ( 2 ) Pv+2kfl{x), X > 0. (1) 

That is, p v \ is a Poisson mixture of central x 2 densities. When v is a positive integer, the non- 
central x 2 arises naturally as the distribution of quadratic forms in normal variables. Specifically, 
suppose Xi, . . . , X v are independent normal random variables with means fik, k = 1, . . . , v, and 
unit variances, then the density of Ylk=i^k * s precisely p U: \ given by (pQ) with A = X^aUi^I- 
The definition ([I]), however, does not require v to be an integer. 

Distributional properties of the noncentral x 2 can be difficult to obtain because the density 
is "not in closed form". Another expression for p u ,\, not necessarily more tractable, is 

PuAx) = ie-(^)/ 2 (|) M/4 V 2)/2 (V^). (2) 
Here I v {x) denotes the modified Bessel function of the first kind given by 

The ratio 

Tu{x) = (3 ) 

plays a crucial role in our analysis. We use I' v (x) and r' v (x) to denote derivatives with respect 
to x. 

We are concerned with shape properties of the noncentral x 2 such as unimodality, log- 
concavity and log- convexity. A function f(x) is unimodal if there exists x* such that f(x) 
increases for x < x* and decreases for x > x* (x* is called a mode). A nonnegative function / is 
log-concave (respectively, log-convex) if log / is concave (respectively, convex). It is well known 
that log-concavity implies unimodality. Since we focus on probability densities supported on 
(0,oo), all statements concerning monotonicity and log-convexity refer to this interval. 



3 Relevant Literature 

Johnson, Kotz and Balakrishnan (1995) discuss the noncentral x 2 distribution theory in great 
detail. Saxena and Alam (1982) study estimation methods for the noncentrality parameter. Ding 
(1992) and Kniisel and Bablok (1996) propose methods to compute the density and distribution 
functions numerically. Concerning distributional properties, Karlin (1968) has proved that if 
v > 2 then the noncentral density p Ut \ is log-concave. Das Gupta and Sarkar (1984) and 
Finner and Roters (1997) derive related results, including log-concavity of the distribution and 
survival functions and generalizations to noncentral F and t distributions; see also van Aubel 
and Gawronski (2003). 

Not only is the noncentral x 2 distribution important for power calculations in hypothesis 
testing, it also has engineering applications including radar detection and communications over 
fading channels. It is intimately related to the generalized Marcum Q function Q u {a, b), which is 
of great interest in information theory and communications. Specifically, Q u (a, b) is the survival 
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Table 1: Values of \ u for selected v G (0, 2). 
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function of a noncentral x 2 distribution with 2u degrees of freedom and noncentrality parameter 
a 2 evaluated at 6 2 : 

/•oo 

Q„(a,b) = / p2 Ua 2(x)dx, a,b>0. 
Jb 2 

See, for example, Marcum (1960), Nuttall (1975), Shnidman (1989), Helstrom (1992), Ross 
(1999), Simon and Alouini (2003), Li and Kam (2006), Sun, Baricz, and Zhou (2010), and 
Yu (2011). Sun et al. (2010) also derive log-concavity results. Yu (2011) establishes the log- 
concavity of Q u (a,b) in b for v > 1/2, confirming a conjecture of Sun et al. (2010). This work 
uses methods that are similar to those of Yu (2011). 



4 Main Results 

For v G (0, 2) define 

A — 2 

9u(X) = r v/2 (y/\(\ + u-4)) - =, \>4-v, 

a/A(A + v - 4) 

where r u (x) is given by ([3J) . We rely on a sign property of g^(A). 

Lemma 1. For each v G (0, 2) there exists a unique \ u G (4 — v, oo) such that g u {\) > or < 
according as A > \ v or < \ u . 

In other words, as A increases, g u (X) crosses zero exactly once from below. The unique 
crossing point X u can be computed by bisection or functional iteration. A small sample is 
recorded in Table 1. Proposition [T] recognizes the limiting cases. 

Proposition 1. We have 

limA I/ = 4, and limA I/ = 2. 

Theorem [H our main result, says that a noncentral x 2 density is unimodal unless the degree 
of freedom is fewer than two and, at the same time, the noncentrality parameter is too large. The 
\ u in Lemma Q] is the critical value that separates the decreasing from the bi- modal categories. 

Theorem 1. The noncentral x 2 density with u > degrees of freedom and noncentrality pa- 
rameter A > is 

i log- concave iff v > 2; 

ii strictly log- convex for small x and strictly log-concave afterwards iff < v < 2 and A > 0; 
Hi decreasing iff < v < 2 and A < \ v ; 

iv bi-modal iff < v < 2 and A > X u . 
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Figure 1: Noncentral x 2 densities with v degrees of freedom and noncentrality parameter A. 

Part (i) of Theorem Q] is well known (see Karlin 1968 or Finner and Roters 1997). But we 
are unaware of any previous work offering the precise distinction between (hi) and (iv). Figure 1 
displays some noncentral x 2 densities, including a bi-modal case (y = 1, A = 5). 

An extreme case of bi-modality is the noncentral x 2 with zero degrees of freedom, which is a 
mixture of a continuous distribution and a point mass at zero. The probability at zero is e~ A / 2 , 
and the density of the continuous part is proportional to ([1]) upon replacing Ylfc=o by SfcLi an d 
letting v = 0. See Siegel (1979) for a detailed exploration of this distribution. 

Next, we study the location of the modes and derive some simple but effective bounds. 
Corollary [T] easily follows from Theorem [T] and its proof. 

Corollary 1. The density p U: \ has 

i a unique mode in (0, oo) if v > 2 or v = 2 and A > 2; 

ii a unique mode at zero if0<u<2 and A < X u ; 

Hi a mode at zero and another in (0, oo) if < v < 2 and A > X u . 

Part (i) of Corollary Q] also appears in Corollary 3.3, part (i) of van Aubel and Gawronski 
(2003). Corollary 1 implies that, if p Ut \ has an interior mode, then so does p u ^ for any A > A. 
Moreover, according to Proposition [21 the mode moves to the right as A increases. This is a 
known result in the log-concave case. 

Proposition 2. // v > 2 or < v < 2 and A > X u then an interior mode M{v, A) exists and, 
for fixed v, M(v, A) increases with A. 
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It is known that M{y, A) increases in u G (2, oo) for fixed A > (van Aubel and Gawronski 
2003). Such monotonicity with respect to v does not seem to extend to the bi-modal case. In 
fact, inspection of Table 1 reveals that, for certain values of A, the interior mode may disappear 
and then reappear as v increases in (0,2). 

When an interior mode exists, Proposition [3] gives more information about its location. 

Proposition 3. // an interior mode M(v, A) exists then it satisfies 

A + I/-4 < M(u,X). (4) 

Moreover, if v > 2 then 

{v-2) M + <M{y,\) < X + u-2. (5) 

If u > 3 then 

A + I/-3 < M(u,\). (6) 

If0<u<2 and A > X u then 

M(y,X) < X + v-3. (7) 

Sen (1989) derives © and @ for v > 2. Van Aubel and Gawronski (2003, Theorem 6.1, 
part (i)) present ([6]) for v > 4. The bounds ^ and ([7]) extend existing results to the bi-modal 
case. The two bounds in ([5]) become tight as A ^ and in fact coincide for A = 0. The bounds 
([6]) and ([7]) are remarkable in view of the asymptotic formula 

M(v,X) = X + v -3 + OiX- 1 ), as A ^ oo. 

We emphasize that M(i/, A) refers to the interior (nonzero) mode. 

In Section 5 we prove the main results. The method is elementary although it does rely on 
properties of modified Bessel functions. Part of the argument may be helpful in establishing 
inequalities in other contexts. Of course, it would be interesting to derive similar results for 
noncentral F or i distributions. 



5 Proofs 

Define l(x) = \ogp Vt \(x). Let us calculate some derivatives in preparation for the main proofs. 
By ([2]) we have 

l' (x ) = -I + v —l + 

2 4x 2 v ^/ ( „_ 2)/2 ( v /A^) 

1 V-2 \[X . r— s 

where (HI) follows from ([3]) and the formula (Abramowitz and Stegun 1972, (9.6.26)) 

4_x(a0 = I v {x) + V -^I v -x(x). 

x 
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We differentiate ([8]) and get 



1 " {X) = + ^V2^) " ^2V2(^) 

2 — v A vyX , r- — . A 2 , /- — . . . 

= + A~x ~ i^V 2 (V^) " ^ /2 (VXx) (9) 

l'(x) l'(x) + 1 - — — (10) 



4x 4 ' \ 2x 

where Q holds by the identity (see, e.g., Amos 1974) 



r' v (x) = 1 - V ' r v {x) - rl(x), (11) 



and (HDD holds by ©. Differentiating ([TDD y ields 

r(x) = - (a + „ _ 4 + 2(i/ - 4)Z'(x)) - /"(x) (V(x) + 1 - . (12) 

We also use the following results, which are easily derived from modified Bessel function 
asymptotics (see, e.g., Abramowitz and Stegun 1972, (9.7.1)). 

x x^ 

r„/2(x) = 27 — r^ +o z;3 )' asx|0; (13) 

' v v z {v + 2) 

r W2( x ) = 1 ~ -7T— 1 + asx-^oo. (14) 

' 2x 

Proof of Lemma[J\ Let v G (0, 2). As A 1 (4 — z/), we have <7„(A) — > — oo by (fT3"P ; as A — >• oo we 
have A^(A) — > 1/2 by {HP - Thus g u {\) changes signs at least once, from negative to positive. 
For any A* such that g v {\*) = 0, we have r u / 2 (t) = i A * ~~ 2)/t, where i = a/ A*(A* + v — 4). A 
calculation using (fTl~D yields 

_ A, + 2v - 6 

r v/2\ t ) ~ 7J 

Thus 

W_ / / 2A, + t/-4 \ 1 / 2A. + 1/-4 

9A " ^ 2(t) V 2^ J " 7 + (A* - 2) ^ ^ 

_ VA. + v - 4 

,3/2 > U - 

A* 

Hence every sign change of g v {\) is from — to +. There cannot be two or more sign changes, 
because one of them has to be from + to — . □ 

Proof of Proposition d We only prove lrnvM A„ = 2; a similar argument works for lim^o A u = 4. 

For fixed A > 2 we have 

A 



hm^(A) = r x (VA(A-2)) - y — . (15) 

By (|13|) . the right hand side of (|15|) is positive for all sufficiently small A > 2. For such A we 
have g v {\) > if v is sufficiently close to 2, and \ u < A by LemmaCQ That is, limsup^ 2 A v — A - 
The claim follows by choosing A I 2. (Obviously \ u > 2 and hence liminf^ \ u > 2.) □ 
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Proof of Theorem [TJ We only prove the "if" parts because the "only if" parts hold after exam- 
ining certain boundary cases and, if necessary, invoking the "if" parts. For example, to show 
the "only if" statement in part (iii), note that if v > 2, or v = 2 and A > X u = 2, then by © 
we have l'(x) > for small x > 0. Hence the density cannot be decreasing. 

Part (i) can be proved by showing that r v / 2 {x)/x decreases in x G (0,oo) (see Saxena and 
Alam 1982 or Yu 2011). By ([8]), if v > 2 then l'{x) decreases in x, as required. 

To prove part (ii) let us assume < v < 2 and A > 0. As x J, we have, by (fl"3j) and ([9]), 
x 2 l"(x) (2- v)/2 > 0. As x -> oo we use (HU) and © to obtain x 3 / 2 Z"(x) -> -VX/4 < 0. 
Thus l"{x) crosses zero at least once from above. If it crosses zero more than once, then there 
exists x* such that Z"(x*) = and l"'(x*) > 0. By (fT2l) we have 

Putting this in (fT0|) we get 



... . A + z^-4 1 ,,. .... , , A + v — 4 (v — 4\ 

4-^^ + 1 ) + ^^(^7j 

+ 5) to- 



phus equalities must hold and we have 

A + ^~ 4 -l'( x )--i 
2(4-!/) 2' 

which yields A = 0, a contradiction. It follows that l"(x) has exactly one sign change on (0, oo), 
and at the change point I'" < 0, yielding a strict (+, — ) sign pattern. 

Let us now assume < v < 2 and A < X u and prove the "if" part of (iii). The argument for 
v = 2 is similar, and the case A = \ v follows by taking a limit. As x 1 we have l'{x) — > — oo; as 
x — > oo we have l'{x) —> —1/2 by ([8]). If does become positive, then there exists < x\ < x 2 
(entry and exit points) such that l'(xi) = l'[x2) = and l"(x\) > > l"(x2)- By (fTUI) we have 
x\ < A + v — 4 = xq < X2, which implies xq > 0. By Lemma [T] 

l'(x ) = ^Lg v {\) < 0. (16) 



Xq 

We can rule out l"{x2) = because it implies x 2 = xq by (fT0|) which contradicts (fT6|) . Because 
/'(x ) < 0, l'(x 2 ) = 0, l"{x 2 ) < 0, there exists x* G (x ,x 2 ) such that /'(x*) = and /"(x*) > 0. 
By (|10p . however, we obtain x* < xo, which is a contradiction. 

Finally, suppose < v < 2 and A > \ u . Phen Z(x) — > oo and Z'(x) — > — oo as x i 0. So one 
mode is at zero. In view of part (ii), let x G (0, oo) be the unique solution of Z"(x) = 0, so that 
l"(x) > or < according as x < x or > x. If l'(x) < then by (fTUI) we have x > A + z/ — 4 = xo- 
(We always have l'(x) + 1 - (u - 4)/(2x) > in ([TO]) by ©.) By Lemma □ 

> l'(x ) = -£=g v {\) > 0, 



which is a contradiction. Phus l'(x) > 0, implying that Z(x) is convex with an interior minimum 
in (0, x) and concave with an interior maximum (the second mode) in (x, oo). Phis proves part 
(iv). □ 
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Proof of Proposition^ Existence of the interior mode follows from Corollary [TJ To show mono- 
tonicity of M(u, A), we prove that h(x) = xr v (x) increases in x 6 (0, oo) for all v > 0. This 
implies that l'{x) as given by (|8j) increases in A, which quickly yields the claim. We have 
h'(x) = r v (x) + xr' u (x). Suppose h'(x*) = for some x* > 0. Then r' u (x*) = —r u (x*)/x* and, 
by QID, rl{ Xif ) + {2v - 2)r u { Xif )/x* - 1 = 0. We get 

^"(x*) = 1 - (2v — 2)r' u (x*) - 2x*rv(x*)r^(x*) - r^(x*) 
= 2 

after simple algebra. Hence any sign change of h'{x) is from negative to positive. However, as 
ijOwe have h'(x)/x — > 1/v > 0. This precludes possible sign changes. Hence h'(x) > 0, as 
required. □ 

Proof of Proposition^ Assume an interior mode exists. We refer to {(z/, A) : v > 2 or v = 
2, A > 2} as "the log-concave case" and {(v, A) : < v < 2, A > A„} as the bi-modal case. 
Denote xq = A + v — 4. If xq < then ([!]) is trivial. If xq > and 1'(xq) > then xq < M(v, A) 
because the density declines after M(u, A). Suppose xo > and 1'(xq) < 0. By (fT0|) . 1"(xq) > 0, 
implying a bi-modal case, and xo belongs to the declining phase in the log-convex region. By 
Theorem [1] we have xq < M(v,\). Thus (J2|) holds. 

It remains to prove ([6]) and ([7]) because Sen (1989) has established ([5]). Denote z = X + u — 3. 
In the case of v > 3 we need to show l'(z) > 0. In the bi-modal case we need l'(z) < instead. 
This suffices because, letting m be the local minimum of Z(x), we have m<xo = A + ^ — 4, 
which follows from l'(m) = 0, l"(m) > and (|10p . Hence m < z, and (J7J) holds if /'(z) < 0. By 
®, 

|/'(,) = V2 (v / Ai)-^i = / l (A). 

Let us follow the proofs of Lemma [Hand Proposition [2] in analyzing the sign pattern of /i(A). 
Suppose < v < 2 and A > 4 — v. As A I (4 — v) we get 

/i(A) -> v^v 7 ?^) - 4r=- ( 17 ) 

V4 — ^ 

One can verify numerically that the right hand side of (I17p is negative for all < ^ < 2 
(there is only one variable over a small range). Suppose there exists some A* G (4 — z^, oo) such 
that /i(A*) = 0. Then r v / 2 {y/\*z*) = (A* - l)/\/A*z* where = A* + z^ - 3, and by {H]), 
^ 2 (\/A*2*) = (f — 2)/(A*z*). Some algebra yields 



r 



2K + V-2, (1 - i/)A* + 3 - z/ 



^.^(V^^^ ^ <18> 
<0. (19) 



2^ 

Thus any sign change of h{\) is from + to — . Since h{\) is negative as A J, (4 — z/), there cannot 
be any sign change. This proves l'{z) < for A > \ v > 4 — v. 

Finally, let us assume v > 3. The calculation (|18p -(|19 p now shows that any sign change of 
/i(A) is from — to +. But h(X) — ^ 00 as A \, 0. Thus sign changes are impossible and we conclude 
that l'(z) > 0, as required. □ 
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Remark. To prove an inequality h(X) > for all A, we show that (i) it holds for small 
(respectively, large) A; (ii) h'(X) is positive (respectively, negative) assuming h(X) = 0. This 
simple method works well for functions involving r u (x) and seems applicable to other problems. 
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